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ABSTRACT 

We present two sets of grid-based hydrodynamical simulations of high- velocity clouds (HVCs) trav- 
eling through the diffuse, hot Galactic halo. These HI clouds have been suggested to provide fuel for 
ongoing star formation in the Galactic disk. The first set of models is best described as a wind-tunnel 
experiment in which the HVC is exposed to a wind of constant density and velocity. In the second 
set of models we follow the trajectory of the HVC on its way through an isothermal hydrostatic halo 
towards the disk. Thus, we cover the two extremes of possible HVC trajectories. The resulting cloud 
morphologies exhibit a pronounced head-tail structure, with a leading dense cold core and a warm 
diffuse tail. Morphologies and velocity differences between head and tail are consistent with observa- 
tions. For typical cloud velocities and halo densities, clouds with HI masses < 10 45 M© will lose their 
HI content within 10 kpc or less. Their remnants may contribute to a population of warm ionized gas 
clouds in the hot coronal gas, and they may eventually be integrated in the warm ionized Galactic 
disk. Some of the (still over-dense, but now slow) material might recool, forming intermediate or low 
velocity clouds close to the Galactic disk. Given our simulation parameters and the limitation set by 
numerical resolution, we argue that the derived disruption distances are strong upper limits. 



Subject headings: Galaxy:halo — Galaxy:evolution 
methods : numerical 



ISM:HI — hydrodynamics — turbulence 



1. MOTIVATION 

Galaxy disks require gas to fuel their ongoing star for- 
mation, yet the fueling mechanism and the nature of 
the gas reservoir remain unclear. In galaxy formation 
models, the halo gas is initially heated as it follows the 
dark matter and subsequently cools and falls to the disk 
(e.g.. iWhite fc ReeslH978l : iWhite fc FrenklH99l . Even 
at z = 0, halo gas is proposed to gradually condense out 
into clouds, and a significant amount (~ 4 x 10 10 M©) 
of hot gas is thought to remain in an extended (~ 
150 kpc), hot, diffuse halo me d ium dMaller fc Bullock! 
2004 iFukugita k Peebles] 120061: iKaufmann et all 120061 : 
Sommer-Larsenl [2006D . Other sources of potential fuel 



for galaxy disks include gas stripped from satellites, and 
possibly some mixed-in galactic fountain gas. Under- 
standing the properties of the diffuse hot halo and the 
clouds that fuel galaxy disks is a crucial step towards 
comprehending how galaxies form and evolve. 

Observationally we know halo gas exists today around 
galaxies th at is likely to fall towards the disk (vg as < 
v^r) (e.g., iFraternali et all 120011: iPutman et all 120021 : 
Thilk er ;t alT 20041 ) . The neutral hydrogen halo clouds 
found around the Milky Way are called high-velocity 
clouds (HVCs) and are presumably a future source of 
star formation fuel for our Galaxy. HVCs are moving 
through a more diffuse halo medium as evident from their 
structure ( e.g.. iBruns et al.ll2000l : IPutman et aLll2003at 
iPeek et all I2007J ) , and from detection s of high- velocity 
OVI absorption (|Sembach et al.|[2003l ). The clouds are 
disrupted by this interaction and form a "head-tail" 
structure, i.e. the head of the cloud is compressed around 
a cold core, and a warmer, diffuse tail extends behind the 
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cloud. The properties and lifetime of a given head-tail 
cloud depend on the density of the diffuse halo medium. 
Head-tail clouds therefore present the opportunity to 
study the lifetimes of halo clouds and the density of the 
diffuse Galactic halo, a medium that is difficult to detect 
directly. 

In this paper we present three-dimensional grid simu- 
lations of HVCs interacting with a diffuse halo medium 
and compare the results to neutral hydrogen observa- 
tions. We discuss two sets of numerical experiments, one 
best described as a wind-tunnel experiment in which the 
model cloud is exposed to a wind of constant density and 
velocity, and a second set following the model cloud on 
its way through an isothermal, hydrostatic halo towards 
the disk. All models include equilibrium heating and 
cooling. The characteristic timescales and pathlengths 
for disruption are assessed in the context of fueling the 
Galactic disk with the observed population of HI HVCs. 
Our main results are summarized in Figure 6, predict- 
ing the distances HVCs of a given mass are expected 
to survive their disruption at a given cloud velocity and 
background halo density. 

2. TECHNICAL ASPECTS 
2.1. The Numerical Scheme 

Calculations were performed with Proteus, an 
unsplit, second-order accurate gas-kinetic scheme 
(iPrendergast fc XdflQQl iSlvz fc Prendergastl H9991 : IXul 
l2001t iHeitsch et al.ll2008[) based on the Bh atnagar- Gross- 
Krook formalism (Bhatnag ar et alJ ll954). The scheme 
conserves mass, momentum and total energy to machine 
accuracy, and it allows the explicit control of viscosity 
and heat cond uction, albeit at a Prandtl number of 1 
(see IXul [2001D . The dissipative constants are set such 
that they c ontrol the (very low) intrinsic numerical dissi- 
pation (see lSlvz et al.ll2002l §3 for a detailed description). 

The thermal physics are implemented as a modification 
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of the total energy via additional heating and cooling 
terms. In the absence of heating and cooling, the gas is 
adiabatic with an exponent of 7 = 5/3. The heating and 
cooling terms are accumulated in a tabulated "cooling 
function" (i.e. energy density change rate e), depending 
on density n and temperature T as e = nT — n 2 A(T). 
The energy change is done iteratively, i.e. if the local 
cooling timescale is smaller than the current Courant- 
Friedrich-Levy (CFL) timestep, we subcycle on the en- 
ergy equation in that cell. This is generally the case in 
the cold dense gas. Heating and cooling are operating 
throughout the whole simulation volume. We derive the 
cooling function from the rates given by IWolfire et al.l 
(|1995af ) for T < 10 4 K, while for higher tem peratures we 
use the rates of Sutherl and fc Dopital (Il993f). We assume 
a metallicity o f 10% solar (jWakker fc van Woerdenlll997t 
lWakkedl200ib . 

The heating is due to a generic UV radiation field, 
which is set t o Go = l.lX(z), w here X(z) = 1/[1 + 
(z/8.53kpc) 2 ] (jWolfire et al.lll995bl ) with z = 10 kpc as 
a "typical" HVC distance. Although we do not know 
the actual distance to each of the head-tail clouds, we 
adopt 10 kpc for most models based on the recent work 
constraining the distance s to many HVCs (|Thom et al.l 
I2006t IWakker et aTll2008ft . Many of the head-tail clouds 
are in the vicinity of larger complexes, suggesting that 
they may have a similar distance. This will be described 
further in the paper presenting the population of HIPASS 
head-tail clouds (Putman et al., in preparation). The 
UV-field is kept constant during the evolution of the 
HVC. All heating and cooling processes are local, and we 
do not address radiative transfer effects. Thus, the clouds 
cannot self-shield, and possible effects of increased UV 
radiation once the cloud approaches the Galactic plane 
cannot be addressed. 

2.2. Setup and Parameters 

We ran two sets of models. Series W are wind-tunnel 
experiments, while in scries H, the cloud falls vertically 
towards the disk's midplane, traversing an isothermal 
halo medium in hydrostatic equilibrium. In both cases, 
the clouds are initially at rest, located in the lower quar- 
ter of the simulation domain. 

All model clouds have an exponential density profile 
(e.g. iBurton et ai1l2001l see Figure []} , convolved by a 
tanh shifted to the nominal cloud radius to make the 
cloud compact. We also experimented with uniform den- 
sity clouds and found that they tend to disrupt even 
faster than the centrally peaked clouds, as expected. 
Wind-tunnel experiments have simulation domain sizes 
varying between 150 2 x 300 and 500 2 x 1000 pc 3 , while 
the free-fall experiments use box sizes of 500 2 x 2000 pc 3 . 

2.2.1. Wind Tunnel Experiments 

For series W, the background medium has a constant 
density and t emperature Tn such that the pressure 
is 100 K cm" 3 (jWolfire et al.lll995ri l for models with 
nh < 10~ 4 cm~ 3 , and 300 K cm~ 3 for models with 
rih = 3 x 10~ 4 cm~ 3 in order to avoid the peak of the 
cooling curve at T m 3 x 10 5 K. The cloud is initially in 
thermal pressure equilibrium with its surroundings. To 
keep the cloud within the simulation domain as long as 
possible, we balance a constant gravitational acceleration 



g = gz against a "wind" entering the bottom of the do- 
main. The wind speed is ramped up linearly at a rate of 
13 km s" 1 Myr -1 until the desired "cloud" velocity vq is 
reached. The velocity vq sets the size of the gravitational 
acceleration g, via 



M c g = TzRlC DPh v 2 J2, 



(1) 



with the cloud mass M c , the drag coefficient Co = 1-0 
and the halo density ph = n^mH. We note that while 
equation (TTJ implicitly assumes that the clouds reach a 
terminal velocity t>o, we use it only as a means to achieve 
an approximate balance between gravity and ram pres- 
sure. It is not clear whether HVCs generally reach their 
term inal velocities ([Benjami n fc DanTvlll997t iPeek et alj 
|2007| ). Moreover, equation (TTJ assumes that the clouds 
do not disintegrate, leading to an over-estimate of the 
terminal velocity vq. The lower z-boundary is defined by 
the "wind" , while all other boundaries are left open. 

2.2.2. Free-fall Experiments 

For series H, we drop the HVC in a hydrostatic isother- 
mal halo at To = 10 6 K. We use a gravitational accel- 
eration of g = 10~ 8 cm 2 s~ 2 , corre sponding to the fit 
betwe en 1 and lOkpc provided by Benjamin fc Danlvl 
( 1997), and we assume a halo density of rih = 10 -4 cm~ 3 
at z = lOkpc. This results in a density profile of 



n h (z) = n exp(-gm H /(k B T )) 



(2) 



with a mid-plane density of no = 4x 10~ 3 cm -3 . The fact 
that this is lower than the volume density of the ionized 
gas in the Galact ic midplane by a factor of 5 — 10 (see 
e.g. lFerrierell2001l ) is of no concern here, since our simple 
model cannot describe the density structure for z < 2 kpc 
correctly anyway. We will discuss in ij4.2l how any sys- 
tematic underestimate of the background halo density 
affects the disruption of the clouds. To keep the cloud 
within the simulation domain, the grid is shifted by ex- 
actly one cell to lower z once the center of the HI mass 
distribution has fallen at least by one cell. The new layer 
of cells at the bottom of the domain is set according to 
the exponential density profile determined by the above 
parameters. This prevents unnecessary and diffusive in- 
terpolation and it allows us to follow the cloud on its trip 
through the halo. 

Apart from this "comoving" grid, the major difference 
to series W is that the halo pressure is no longer constant. 
This leads to a successive compression of the cloud and 
its fragments, until they drop below the resolution limit 
of th e dynamical instabilities r elevant for cloud disrup - 
tion (|Roulston fc Ahrend[i997l : iKorvcanskv et al.ll2000h . 
This effect is amplified by the increasingly higher cooling 
rates. In the extreme case, the remaining cloud fragment 
can collapse to a thin "needle" in essential free-fall. To 
prevent such unphysical behavior, we regrid the simula- 
tion whenever the most massive fragment drops below 
the resolution limit. For such models, we simply reduce 
the domain size by one half, at the same number of cells, 
thus doubling the spatial resolution. Typically, one re- 
gridding step is necessary to follow the cloud evolution 
until disruption. It should be noted that even with the 
regridding the smallest, cold dense fragments are not re- 
solved with respect to the hydrodynamical instabilities. 
Thus, their lifetimes tend to be over-estimated. 
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is traveling nearly head-on towards the observer. Model 
Wblal5b at 45° develops a pronounced head-tail struc- 
ture reminiscent of the structures in the observed cloud, 
including several fragments in the tail. Shear flow insta- 
bilities act along the cloud and tail flanks, and the tail 
itself shows kinks as a result of a turbulent wake. Cool- 
ing leads to fragmentation, resulting in "multiple cores" 
along the tail. In the head-on view (cloud Wblal5b, 
5°), the tail is nearly completely suppressed, but might 
be discernible in a velocity shift between the cold and 
warm gas component. We will explore the kinematic 
signatures of the clouds and their disruption in a subse- 
quent paper. The panel on the right (model Welc05b, 
30°) shows a more massive cloud at a late evolutionary 
stage. The cloud has disintegrated into multiple cores. 
The minimum temperatures in the cores reach « 100 K. 



Fig. 1. — Logarithmic density profiles for models as indicated. 
The profiles are set by an exponential convolved with a tanh shifted 
to the nominal cloud radius to make the cloud compact. 

2.2.3. Model Names 

The model parameters are summarized in TableQ] The 
model names consist of seven characters, the first indicat- 
ing whether the model is a wind-tunnel experiment (W, 
ffT2~Tj) or a free-fall experiment (H, §2~2l>]) . The sec- 
ond character denotes the cloud radius (a = 21, b = 50, 
c = 89, d = 133, e = 150 pc). The initial peak cloud den- 
sity is given by the third character in 0.1 cm -3 , while the 
fourth character indicates the halo density (o = 3x 10~ 5 , 
b = 1 x 10~ 4 , c = 3 x 10~ 4 cm" 3 ). The fifth and sixth 
character stand for the nominal cloud velocity for wind- 
tunnel experiments, and for the peak cloud velocity for 
free-fall experiments, in 10 km s . Finally, the seventh 
character stands for the resolution along the shorter axes 
(a = 128, 6 = 256 cells). 

3. RESULTS 

We start our discussion of the results in terms of 
morphology fi )3.ip . including a comparison to observed 
HVCs. A detailed view of three representative models 
is given in §3.21 We also summarize the cloud life times 
and their travel distances ( N3.3p . discuss resolution ef- 
fects f £|3.4[) , and differences to previous work ( £|3.5i see 
also g2Z3). 

3.1. Morphologies 

Figure [2] shows column density maps of HI (i.e. gas at 
T < 10 4 K for the model clouds) of observed and modeled 
HVCs with the same column density contours for each. 
Observed clouds were taken from t he HI Parkes All-Sk y 
Survey (HIPASS) cloud catalogue ()Putman et alJl2002f ). 
Clouds with a head-tail structure were found for ~ 40% 
of the small HVCs (CHVCs and :HVCs) in this catalogue 
(Putma n et al., in prep.), and a lso along the Mage llanic 
stream (|Putman et alj l2003bf ) . iBriins et all (j2000f ) put 
together a flux-limited sample of head-tail clouds with 
the Leiden-Dwingeloo Survey (35' resolution vs. 15.5' for 
HIPASS) and found 20% of their entire sample show a 
velocity and column density gradient consistent with be- 
ing head-tail clouds. The modeled clouds in Figure [2] 
have been degraded to a "beam width" of 3.0', and they 
have been rotated by the amount indicated on top of 
each model map. For small rotation angles, the cloud 



3.2. Cloud Evolution and Disruption Mechanisms 

We will discuss the cloud evolution in terms of 
three models, namely Wblal5b, Hclbl3a and (at 
twice the resolution, but otherwise identical parameters) 
Hclbl5b, representing windtunnel and free-fall models, 
and demonstrating resolution effects as well as long- 
term evolution. Figure [3] summarizes the time history 
of the three models, from top to bottom the mass, veloc- 
ity, Mach number, ratio of thermal over dynamical time 
scales and the number of grid cells against time. We 
will discuss each of them in turn. To measure the time 
histories of the clouds, we ran a core finder identifying co- 
herent structures in three dimensions with temperatures 
T < 10 4 K as a proxy for neutral hydrogen. Structures 
with fewer than 64 cells are not counted. In the follow- 
ing, we will use "cloud mass" synonymously with "cloud 
mass in HI" . 

3.2.1. Mass evolution 

We begin our discussion with the mass histories (top 
panels of Fig. [3} , which indicate the survival chances of 
the model cloud. The simulations start out with one cold 
cloud, which subsequently fragments due to a combina- 
tion of dynamical and thermal instabilities. The thick 
solid line denotes the total mass in HI, while the thin 
solid line follows the mass evolution of the most massive 
fragment (usually the main cloud except for very late 
stages of nearly complete fragmentation). The dashed 
line indicates 10% of the initial cloud mass. The center- 
of-mass velocity of each fragment is shown in color. 

With increasing time, more and more material is ab- 
lated from the cloud, resulting in a growing number of 
cloud fragments. The increasing lower mass envelope of 
the fragments in model Hclbl5b (to some extent visible 
in Hclbl3a as well) results from the recooling of the frag- 
ments. With decreasing z, the cloud fragments are being 
compressed by the increasing background pressure, thus 
leading to higher cooling rates and "re-formation" of cold 
Hl-cloudlets. The mass fraction in the small, "re-formed" 
cloudlets amounts to w 10% of the Hl-mass of the origi- 
nal cloud. This effect is absent in the wind-tunnel exper- 
iments, where the cloud is being shredded into succes- 
sively smaller fragments. The evolution of "real" HVCs 
we expect to lie in between these two extremes of expo- 
nential density increase and constant density of the back- 
ground medium, since the cloud trajectories will gener- 
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Cb3c, 45deg Cb3c, 5deg CbOa 3Udeg 




Fig. 2. — Column density maps of observed and modeled HVCs. Observed clouds were taken from the HIPASS catalogue by Pu tman et aEI 
(2002). From left to right: HIPASS cloud #648, model cloud Wblal5b at 45 Myr, rotated by 45°, HIPASS cloud #200, model cloud 
Wblal5b at 45 Myr, rotated by 5°, and model cloud Welc05b at 110 Myr, rotated by 30°. Contours indicate column densities levels 
identical for all clouds, at (1.2,1.6,2,4,8,12,16,20) X 10 18 cm -2 . Because of the higher peak densities, model Welc05b has additional 
contour levels at (40,80, 160) X 10 18 cm -2 . The models have a nominal resolution of fs 2 pc, corresponding to 0.7' at 10 kpc. The maps 
have been degraded by a Gaussian corresponding to an angular resolution of 3.0' for a more realistic comparison to observations. 




Fig. 3. — Time histories for models Wblal5b, Hclbl3a and (at twice the resolution) Hclbl5b. From top to bottom: logarithm of Hi 
mass of all cloud fragments, their velocity, the internal Mach number of each cloud fragment, the ratio of cooling time scale over dynamical 
timescale to estimate the likelihood for disruption, and the number of cells in each fragment. Color codes as indicated in the panels. Note 
that the velocities are given in 10 km s~ 1 . The ratio Tc/t^ is color-coded with the temperature as in the Mach number plot. Solid lines 
denote the evolution of the most massive fragment (i.e. the main body of the cloud). For the mass-history, a thick solid line shows the 
total Hi mass. The dashed line at the bottom of the cell number plot indicates the threshold of 64 particl es re quired to be accepted as a 
"cloud". Note that this is far below the resolution limit for hydrodynamical instabilities (see discussion in i|3,4l l. 
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TABLE 1 
Model Parameters 



Name N x n c vo Mq Nq Ro Rn to -Do 







cm 3 


cm 3 


km s 1 


1O 3 M 


10 19 cm- 2 


pc 


cells 


Myr 


kpc 


WalblOb 


256 


i x icr 4 


0.1 


100 


1.4 


1.32 


21 


36 


35 


3.3 


Wa3bl0b 


256 


l x icr 4 


0.3 


100 


4.5 


3.96 


21 


36 


65 


5.9 


Wblal5a 


128 


3 x 10~ 5 


0.1 


150 


5.9 


1.32 


50 


32 


50 


5.4 


Wblal5b 


256 


3 x icr 5 


0.1 


150 


3.0 


1.32 


50 


51 


59 


7.2 


Welc05a 


128 


3 x 10~ 4 


0.1 


50 


41 


2.94 


150 


38 


150 


5.3 


Welc05b 


256 


3 x icr 4 


0.1 


50 


41 


2.94 


150 


77 


162 


6.0 


Welc09a 


128 


3 x icr 4 


0.1 


90 


41 


2.94 


150 


38 


81 


5.4 


WelblOa 


128 


1 x 10~ 4 


0.1 


100 


41 


2.94 


150 


38 


64 


3.2 


Hclbl3a 


128 a 


i x icr 4 


0.1 


130 


20 


2.55 


89 


16 


92 


9.0 


Hclbl5b 


256 


1 x 10~ 4 


0.1 


150 


20 


2.55 


89 


92 


68 


9.2 


Hc2bl5a 


128 a 


l x icr 4 


0.2 


150 


40 


5.10 


89 


46 


39 


6.7 


Halb06a 


128 


l x icr 4 


0.2 


60 


7.3 


1.32 


21 


18 


36 


1.8 


Hdlbl7a 


128 a 


1 x 10~ 4 


0.1 


170 


37 


2.87 


133 


34 


103 


11 


Hdlbl6b 


256" 


l x icr 4 


0.1 


160 


37 


2.87 


133 


68 


74 


11 


Hdlb21a 


128 a 


1 X 10~ 4 


0.1 


210 


37 


2.87 


133 


34 


103 


12 



Note. — 1st column: model name (see text), 2nd: number of grid cells along shorter domain 
axes (a = 128, b = 256 cells), 3rd: background halo density (a = 3 X 10 -5 , b = 1 X 10 -4 , 
c = 3 X 10 -4 cm -3 ), 4th: cloud peak density, 5th: nominal cloud velocity, 6th: initial cloud mass, 
7th: initial peak column density, 8th: initial cloud radius (a = 21, b = 50, c = 89, d = 133, 
e = 150 pc), 9th: resolution elements per cloud radius, 10th: characteristic life time (eq. [3]), 
11th: characteristic distance traveled (eq. [3]). 
a These models have been regridded (see i|2. 2.11 1. 
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ally not be restricted to motions along z, specifically if 
the clouds stem from material stripped off satellites. 

3.2.2. Peak Velocities 

The second row of panels in Figure [3] complements the 
mass histories discussed above. Although this informa- 
tion may seem redundant, it is easier to glean the peak 
velocities reached, and the velocity distribution of the 
fragments. As in the mass histories, the thin solid line 
traces the evolution of the most massive fragment. 

The velocities of the cloud fragments are generally 
smaller than that of the main body, indicating a growing 
"tail" of the cloud, due to dynamical drag. The lag is 
consistent with the veloc ity grad ients in observed HVC s 
(|Bruns et al.ll2000Ll200H see also IWestmeier et al.l[2005h . 

While for the wind-tunnel experiment, the velocity is 
prescribed, the velocities reached by the H-series range 
well b elow ballistic velocities (see e.g. [Benjamin" fc Danlvl 
1997), and rarely exceed « 150 km s _1 . Several effects 
play a role here. First, because the clouds are deformed 
due to hydrodynamical interaction with the background 
gas, the drag forces might be higher than expected. Sec- 
ond, the clouds lose mass and eventually disintegrate, 
and third, the increasing background density leads to 
stronger drag. In fact, for model Hclbl3a, the combina- 
tion of mass loss and background density increase brings 
the cloud fragments to a halt at around 100 Myr. To- 
gether with the recooling of cloudlets (see previous sec- 
tion), this offers an interesting possibility of reforming 
Hl-clouds close to the disk 

Observationally, relative to the Galactic standard 
of rest (GSR), the maximum HVC velocities are 
±250 km s _1 , with the more typical GSR velocities 
< 100 km s _1 in ma gnitude and a median value 
of -30 km s " 1 (e .g. IWakker fc van Woerder] 119911 : 
iPutman et aTll2002j ). Though we do not know the true 
velocities of the clouds, this represents our only obser- 
vational constraint. Many people have begun to use the 
deviation velocity, or how much the HVC motions dif- 
fer from a simple model of Galac tic rotation, when dis- 
cussing HVCs (e.g. IWakkerll2004l) . The majority of the 
deviation velocities are < 200 km s _1 . The peak val- 
ues of ±300 km s _1 belong to HVCs associated with the 
Magellanic System. 

3.2.3. Internal Mach Numbers 

The third row of Figure [3] shows the logarithm of 
the internal Mach number for all identified fragments. 
We define the Mach number as the density-weighted 
rms velocity within a fragment over its average sound 
speed. As an estimate of the sound speed, we show 
the volume-averaged temperature of each fragment in 
color. Note that the lowest temperatures can reach 
100K in our models, in a greement with observed values 
(ISchwarz fc Wakkenl200i) . 

The Mach numbers reach values up to 2, although 
it should be noted that these are averaged over the 
cloudlets - in the coldest regions, the Mach numbers can 
exceed values of 5. Thus, the clouds have strong inter- 
nal dynamics, indicating dynamical instabilities at work. 
Yet the Mach numbers and disruption timescales do not 
correlate - for that our measures at this point are too 
crude, and we defer a more detailed discussion to a fu- 
ture paper. 



3.2.4. Thermal over Dynamical Timescales 

The ratio of the cooling time over the dynamical 
timcscale t c /t,j can be used as a measure for the im- 
portance of thermal over dynamical instabilities. Cool- 
ing can stabilize the cloud against disruption by damping 
the p ropagation of waves into the cloud (e.g. iVietri et al.1 
I1997T) . The cooling timescale is given by r c = ksT/ (nA), 
and as the dynamical timescale we take the sound cross- 
ing time D/c s , where the diameter D is estimated by 
taking the geometric mean of a cloud fragment's extent 
along the three grid axes, i.e. D = (DxDyDz) 1 / 3 . For 
Tc/Td > 1, the cooling time is longer than the sound 
crossing time, and perturbations will only be slightly 
damped when traveling into the cloud. For T c /rd < 1, 
thermal effects dominate, leading to a stabilization of the 
cloud. If the densities and temperatures of the cloud frag- 
ments were constant with time, dynamical effects should 
eventually win over thermal ones with progressive frag- 
mentation of the cloud. 

Cloud Wblal5b is completely dominated by dynamics 
{Tc/Td > 1 for all t). This is consistent with the gradual 
disruption of the cloud mirrored by the continuous mass 
loss (top panel). As the number of cells in the most 
massive fragment (bottom panel) stays nearly constant, 
the mean cloud density decreases, thus further lowering 
the importance of thermal effects. 

Clouds Hclbl3a and Hclbl5b arc initially dominated 
by dynamics, however, the increasing background pres- 
sure compresses the cloud and thus shortens the cool- 
ing time, rendering some of the fragments and eventu- 
ally the main body of the cloud thermally dominated for 
t > 45 Myr. Yet the cloud still disrupts because the dense 
gas enters the thermally unstable regime of the cooling 
curve. Fragmentation is enhanced since the condensation 
mode of the thermal instability has an outer length scale 
set by the pr oduct of the sound crossing time and the 
coolin g time (|Burkert fc Lin|[2000t iHeitsch fc Hartmannl 
l2008h . 

3.2.5. Cloud Resolution 

The bottom panels of Figure[3]show the number of cells 
in each cloud (fragment). The dashed line denotes a limit 
of 64 cells below which a fragment is not counted as a 
cloudlet any more. Clearly, this number is well below 
the hydrodynamical resolution limit (see discussion in 
t j3.4[) . For the vast majority of the time, the main body 
of the clouds stays above 10 4 cells. The jump at 60 Myr 
in model Hclbl3a is due to the regridding of the model 
(see 32~2~T|) . 

3.2.6. Cloud Travel Distances 

Figure [5] shows the mass and velocity history in terms 
of the cloud's position z above the Galactic plane, for 
models Wblal5b, Hclbl3a and Hclbl5b. Strictly, for 
model Wblal5b, the initial position is arbitrary. Also 
note that the cloud velocities for the W-series are not 
self-consistent, but are imposed by the wind entering the 
simulation domain (see M2.2.1|) . For the H-series, the ve- 
locity profiles of the most massive fragment (i.e. the main 
bo dy of the cloud) are initial ly similar to those derived 
by iBeniamin fc Danlvl (|1997t ) (their Fig. 3) for clouds 
with a column density of 10 19 cm -2 being dropped at 
10 kpc. Yet our clouds do not reach terminal velocities, 
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but are successively decelerated due to continuous mass 
loss. As demonstrated by models Hclbl3a and Hclbl5b, 
the more massive fragments are to be found at lower z, 
and - as discussed above - the fragments lag behind the 
main body of the cloud. This should be observable as a 
gradient along head-tail clouds or a sh ift in velocity com- 
pone nts between warm and cold gas (|Bruns "etatll2000L 
l200lh . 

3.3. Characteristic Disruption Times and Distances 

Our model data allow us to estimate characteristic 
timescales tq and distances Dq within which a cloud of a 
given mass and at a given background halo density will 
lose its HI content. This information we then can use 
to predict whether observed HVCs with distance con- 
straints will retain any of their HI, or whether they will 
disintegrate before they reach the disk. The time series 
of our models provide us with an opportunity to sample 
a whole range in cloud mass and (for series H) in halo 
densities. 

The disruption time and distance are given by 

T o = ttftii ( 3 ) 



and 



Do 



dM/dt 



M 



(4) 



(5) 



dM/dz ■ 

with the cloud mass M, the mass loss rate with respect 
to time dM/dt and distance dM/dz. The mass range is 
sampled by taking the differences (here for time) 

dM = M(t) - M end 

dt t t en( i 

for all model times t > tend, where t en d is the smaller of 
the time at which the cloud has lost 90% of its original 
HI content, and the end of the simulation. Figure[5]sum- 
marizes Do and tq in dependence of the cloud mass, for 
all models. Colors indicate the velocity v(t) of the cloud. 

Most clouds are disrupted within 10 kpc and 100 Myr, 
with the majority surviving for not longer than 8 kpc. 
Slower clouds survive for longer times and distances, mir- 
roring the lower mass-loss rates, with clouds at v z > 
100 km s -1 lasting for at most 8 kpc. Higher background 
halo densities also lead to faster disruption, as expected. 

Figure [6] can be used to predict the disruption dis- 
tances Do for observed HVCs, given their mass, the halo 
background density rih and the cloud velocity v. The 
latter two have been combined into the ram pressure (in 
K cm -3 ). The left panel of the Figure shows the dis- 
ruption distance D at a given mass and (in grayscale) 
ram pressure of all models, derived from Figure [5] Over- 
plotted are the masses and distances above the plane z 
for a sample of HVCs taken from the extended catalogue 
bv lWakker fc van Woerder] (|199lD with ~ 30' resolution, 
and from a catalog of the clouds at the tail of Complex C 
as observed with Arecibo by Hsu et al. (2009; 3.5' reso- 
lution). We selected only clouds of HVC complexes with 
known distance constraints (see Tab [2|). Masses have 
been calculated from the distance constraints and the 
total flux per cloud as Mjj/[M q ] 

[Jy km s -1 ]. 

Clearly, model clouds with higher mass tend to travel 
farther distances until their HI content is completely gone 



2.36 x 



TABLE 2 

Distance Constraints 



Complex 


D [kpc] 


z [kpc] 


reference 


ACHV/Cohen Strm 


5.0-11.7 


0.9-9.4 


(1) 


GCP/Smith Cld 


9.8-15.1 


2.6-10.1 


(1) 


A 


4-10 


2.2-7.7 


(2) 


M 


1.5-4.4 


1.1-4.0 


(3) 


C 


10.0 ±2.5 


0.9-11.0 


(4,5) 


WB 




1.7-10.2 


(6) 



Note. — 1st column: complex name, 2nd column: he- 
liocentric distance, 3rd column: distance towards Galac- 
tic pl ane, 4th column: reference: (1) IWakkej^etjil 
(2008), (21 Ivan Woerden et aO lTl999T). (3) IDanlv et al 
119931. f4llThom et all (120081 1. (51 IWakker et all (120071 1 



1 
1- 



(6) Thorn et al. (2006). The ranges in z have been es- 
timated using fiducial distances D of 8.4, 12.5, 7.0, and 
3.0 kpc for complexes with D-ranges given, and then 
accounting for the full extent in Galactic latitude. 

- indicated by the upper envelope of the grayscale map. 
This effect is offset by an increased ram pressure - an 
increase by a factor of 10 reduces the disruption length 
scale Do by approximately the same factor. 

Figure [6] shows that many of the observed HVCs will 
not make it to the disk in cold form. The left panel is re- 
stricted to the mass range sampled by our model HVCs. 
Clouds located above the gray shaded region will clearly 
disrupt before making it to the disk, while clouds within 
the Gray shaded region may make it depending on the 
exact halo density they are moving through and their ve- 
locity. Yet since more than 90% of the mass in the HVC 
complexes with distance constr aints resides in clouds of 
large masses (as catalogued by [Wakker fc van Woerdenl 
1 1 9 9 If ) . the right panel of Figure [6] shows an extrapolation 
of the ram pressure map to larger masses. The extrapola- 
tion has been calculated by a two-dimensional quadratic 
fit to the ram pressure surface shown in the left panel. Al- 
though this approach is rather crude, it demonstrates the 
salient correlation between Do, M and P ra m- While the 
smaller HVCs will most likely be disrupted on their way 
to the disk, the situation is less clear for the main body 
of the larger complexes at these distances. Yet it should 
be noted that when higher resolution catalogs are used, 
these massive clouds are o ften bro ken up into smaller 
cloud s (|Putman et al.ll2002t Fig. 1 of lSchwarz fc Wakkeri 
2004). In other words, the fact that the substructure 
in HVC complexes tends to be underestimated renders 
the predicted disruption distances of Figure [6] as upper 
limits. 

Figure [7] offers a straight-forward way 3 to read off the 
disruption length scale Do as a function of log M and the 
ram pressure log P. This has been derived by inverting 
the two-dimensional quadratic fit, yielding 



Do = 3.481 log M 



165.2 + 77.16 log M 

18.54(logA/) 2 + 86.89 log P 
[kpc], 



1/2 



(6) 



where M is given in M©, and P in K cm -3 . Since the 
quadratic form represents the data reasonably well only 
in a limited range of parameters, equation ([6]) is limited 

3 We thank the referee for the suggestion. 



Fate of High- Velocity Clouds 




Fig. 4. — Mass (top) and velocity (bottom) of cloud fragments against distance from Galactic plane for models Wblal5b, Hclbl3a and 
Hclbl5b. Note that the velocities are given in 10 km s~ 1 . Solid lines denote the evolution of the most massive fragment (i.e. the main 
body of the cloud). For the mass-history, a thick solid line shows the total Hi-mass. Because of the small mass contained in the fragments 
of model Wblal5b, the lines are indistinguishable for that model. 



to the range in M and P indicated by the outer con- 
tours in Figure [Jj Within that range, equation (|6|) is 
accurate to approximately 10%. Below the thick line to- 
wards high masses and low pressures (for Dq exceeding 
16 kpc), equation ([6]) breaks down. Still, it allows us to 
read off the approximate disruption length scale Dq for 
physically reasonable values of M and P. 

3.4. Resolution issues 

Numerical resolution is a main concern when mod- 
eling the disruption of clouds by dynamical instabili- 
ties. In their two-dim ensional stu d y of cloud disrup- 
tion by passing shocks, iKlein et al.l (| 19941 ) quote a res- 
olution of approximately 100 cells per cloud radius (or 
i?100) to fully resolve the f low d ynamics around the 
cloud. Yet iKorycanskv et all (|2000t ) noted that since the 
dynamical instabilities lead to non-linearities and tur- 
bule nce, detailed convergence might be unfeasible (see 
also lMac Low fe Zahnlelll994h . They argue for statisti- 
cal convergence (or converge nce of integrated qua n tities) 
at ks i?100. In a similar vein. iRoulston fe Ahrensl (|1997f ) 
argued in the context of the disruption of comets in plan- 
etary atmospheres that i?25 is sufficient to converge on 
the penetration depth within a scale height. 

Checking the lower panels of Figure [3l it is obvious 
that our models will have a hard time meeting the reso- 
lution criterion of i?100 for all of the cloud fragments, 
especially at later times when the main body of the 
cloud will develop a strong tail, reducing its cross section 
and thus possibly suppressing Rayleigh- Taylor instabil i- 
ties (IRoulston fc Ahrensl [l997t IKorycanskv et al.ll2000fl . 
Clearly, the cloud fragments will be dynamically under- 
resolved, an effect which is exacerbated in the models 
of the H-series, due to the increasing external pressure. 
The small est cloudlets might be do minated by evapora- 
tion (e.g. Stanimir ovic et aT1l2008l in the context of the 
Magellanic Stream). 

Figure [5] helps us to assess the effect of Rn on the 
mass loss rate in terms of time and distance traveled 
(eqs. |3I4| ). The mass loss rates will be correlated with 
the strengt h of the dynamical instabilities disrup ting the 
cloud (e.g. lNulsedll98l IKorycanskv et al.ll2000D . which 



in turn will be suppressed in underresolved models. No 
correlation between mass loss rates and resolution is dis- 
cernible, if we ignore the clearly underresolved model 
Halb06a. Although this clearly does not mean that the 
models are fully resolved, we are confident that they are 
sufficiently resolved for our purposes. 

3.5. Differences to Previous Numerical Work 

The wind-tunnel experiments (series W) are conceptu- 
ally si milar to the models discussed by iQuilis fc Moord 
(|2001f h Consistent with their claim, we do not no- 
tice a substantial difference in disruption timescales be- 
tween models with and those without cooling 4 . How- 
ever, for series H, the comparison models without cool- 
ing showed substantially shorter disruption times and 
distances, since the cloud suffers adiabatic compression 
and heating due to the increasing background pressure. 
Thus, the initially cold material is rapidly heated, and 
thus "lost" (i.e. not identified as HI any more). In that 
sense, cooling is relevant for the longer survival of the 
clouds of series H, and it is obviously indispensable for 
t he recooling obs e rved i n e.g. model Hclbl3a. 

ISantillan etaTI (|2004l ) studied the impact of HVCs 
on the Galactic disk using three-dimensional isothermal 
magneto-hydrodynamical (MHD) models. Focusing on 
the physics of the impact, they had the HVC start at 
z = 2.5 kpc above the plane. Thus, conclusions on the 
possibly stabilizing effects of magnetic fields on HVCs 
over long distances could not be addressed. 

A related proble m albeit in a different ph ysical regime 
was considered by [Ruszkow ski et all (|2007l ) who studied 
the disruption of buoyant bubbles driven by AGN feed- 
back in galaxy clusters. They found that tangled mag- 
netic fields can prevent the disruption of the bubbles if 
the coherence length of the field is larger than the bubble 
size, i.e. if the field is only weakly tangled with respect 
to the bubble scale. Otherwise the bubbles are quickly 
s hredded. 

iBland-Hawthorn et al.l (|2007f ) focus on the origin of 

4 Although we ran a few comparison models without cooling, we 
decided not to discuss them in more detail, since their restricted 
physics render them irrelevant for our purposes. 
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Fig. 5. — Distance Do and timescale to within which a cloud of 
given mass will lose its total Hi content. Each "track" represents 
the time sequence of one model, with each point indicating Do 
or ro of the largest fragment (see eqs. 13141 and text). Thus, the 
disruption rate dM/dt increases over time. Colors denote the cloud 
velocities, showing that faster clouds are shorter-lived. Symbol 
sizes indicate model times, with large symbols for early times, and 
small ones for late times. Note that Do is not a distance above the 
plane, but a disruption length scale. 

the Ha emission along the Magellanic Stream and ex- 
plain it by a model of successively shocked material as 
it is ablated from clouds in the Stream predom inantly 
by K elvin-Helmholtz instabilities (see however iMuralil 
2000). Since they are interested in the dynamics of the 
Stream's gas, they start with a two-phase medium, with 
the cold phase having a fractal distribution. Their nu- 
merical method includes non-equilibrium cooling. They 
speculate that much of the gas accretion necessary to 
fuel Galactic star formation happens in the warm ion- 
ized phase, largel y of material shorn off the Magellanic 
Stream (see also iBland-Hawthornl [2008) . Clearly, as 
shown above, HI clouds originating at the distance of 
the Magellanic Stream will not survive the trip to the 
disk. 

4. DISCUSSION 
4.1. Interpretation of Morphologies 



The model clouds usually develop a pronounced tail, 
visible not only in the maps (Fig. [2]), but also in the veloc- 
ity lags of the smaller HI fragments (Figs. 1314ft . The lag 
is consistent with the velocity gradients and/or shifts be- 
tween warm and cold ga s compo nents in observed HVCs 
(|Bruns et alJl2000L l2001t see also IWestmeier et al.ll2005t l 
and should be examined further in higher resolution sur- 
veys with Parkes and Arecibo. The head-tail structure 
will be most pronounced when the trajectory of the cloud 
is perpendicular to the line-of-sight. HVCs that appear 
round may show a kinematic offset in their cold and warm 
components providing additional evidence of clouds be- 
ing destroyed by the surrounding halo medium. In the 
future it should be possible to create a grid of the typi- 
cal length/ velocity offset expected for clouds at different 
stages of destruction and viewing angles. 

4.2. Interpretation of Timescales and Distances 

From Figure [5] we see that the majority of our model 
clouds will be disrupted within < 8-10 kpc. However, 
various physical and numerical issues will affect our esti- 
mates. 

4.2.1. Effects Reducing to and Do 

The following effects will render our estimates as upper 
limits. 

(1) — The actual halo density may be l arger than our 
maxim um of nh = 3 x 10~ 4 cm -3 . iGaensler et al.l 
(2008) set limits on the halo density at z > 5 kpc 
using the 6 measurements they have that fit that cri- 
terion. They find a 3er upper limit of rih < 7.6 x 
10~ 4 cm~ 3 . This type of limit is consistent with 
other work based on pulsar dispersion measures (e.g. 
iHowk et al.l l2006h . The majority of the constraints on 
the halo density are integrated with assumptions on 
the temperature and structure of the halo. To ex- 
lain observations of OVI absorption, ISembach et al.l 
pro p ose d ensities < 10~ 4 down to 10~ 5 cm~ 3 . 



plain 



iPeek et all (|2007h find consistency with a density of a 
few 10 ~ 4 cm -3 at about 10 kpc, and lGrcevich &: Putman! 
(|2009h estimate that a halo density > 10~ 4 cm -3 
is likely at these typ es of distances to explain dwar f 
sph eroidal stripping. [Breeman fc Llovd-Daviesl (p007h 
and lSanders et all \20m combine X-ray absorption and 
emission measures and find a densi ty of 9 x 10~ 4 cm " 3 
for a lengthscale of 20 kpc. See also (Bregmanl (|2007f ) for 
a discussion. 

(2) — The inflow gas density in the W-series is constant, 
i.e. we do not take into account the increase of the gas 
density along the trajectory of the cloud towards smaller 
z for those models. 

(3) — The velocity of our clouds gradually increases to 
the velocity listed in Table [TJ so any initial velocity from 
the origin of the cloud - e.g. the motion of the Milky 
Way satellite - is not included. 

(4) — Resolution effects - in the H-series exacerbated by 
the increased cooling due to the rising external pressure - 
tend to work in favor of faster disrup tion: under-resolved 
clouds tend to take longer to d isrupt feoulston fc Ahrensl 
U997tlKorvcanskv et al.ll2000t see also Fig. |SJ). 




Fig. 6. — Left: Disruption distance Do derived from models shown in Fig. \E\ against cloud mass M, for a range of background halo 
densities and cloud velocities, combined into a ram pressure mjjn^v 2 /kg (see colorbar). A ram pressure of 10 4 K cm -3 corresponds to 
approximately 100 km s _1 at a background halo density of 10 — 2 cm -3 . Symbols represent the z height above the plane and mass of clouds 
that are part of HVC complexes with known distance constraints (see Tab. [2}- Right: Similar to left panel, but using a two-dimensional 
quadratic fit to the ram pressure surface, and extrapolating to higher masses to estimate the disruption distances for more massive clouds. 
Except for the main body of Complex C (at fs 3 X 10 6 Mq), all clouds of the observational comparison sample are shown. 




Fig. 7. — Contours of the disruption length scale Do in kpc, as 
a function of the cloud mass M and the ram pressure P. A ram 
pressure of 10 4 K cm~ 3 corresponds to approximately 100 km s — 1 
at a background halo density of 10 -2 cm -3 . The range of appli- 
cability of equation JBJ is indicated by the 0-line towards higher 
P and lower M, and the thick line towards lower P and higher M 
(i.e. D < 16 kpc). 



(5) — Assuming that the cooling timescale does not sub- 
stantially decrease during the lifetime of the cloud, break- 
ing the cloud into successively smaller fragments should 
render t c /t,i > 1 eventually, i.e. the fragments should be 
more likely to disrupt due to dyn amical instabilities, an d 
eventually through evaporation (jCowie fc McKedll977h . 

(6) — The substructure in HVCs and comple xes tends to 
be underestimated (see iPutman et all l2002t and Fig. 1 
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Fig. 8. — Mass loss rate in terms of time (bottom) and distance 
traveled (top), against the number of resolution elements per radius 
Rn. Ignoring the clearly underresolved model Halb06a, there is no 
trend of mass loss rates with resolution. 



of ISchwarz &: Wakkerl [2004). Existing substructure in 
HVCs (whether imprinted during their formation or 
already due to the interaction with the background 
medium) will accelerate the disruption. The presence 
of substructure is equivalent to a gas volume filling fac- 
tor fv < 1, i.e. some part of the cloud volume is at 
densities lower than the mean density (and if fv <C 1, 
a large part). Since the dynamical instabilities respon- 
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sible for fragmentation have shorter growth timescales 
at lower density contrasts, a small volume filling factor 
promotes cloud disruption. Moreover, the sound speed 
will be larger and the cooling time will be longer in a 
large part of the cloud volume, resulting in a dominance 
of dynamical over thermal effects. 

4.2.2. Effects Increasing tq and Dq 

There are also mechanisms that will potentially in- 
crease the lifetimes of the halo clouds. Magnetic fields 
can s uppress the dynamical instabilities (|Chandrasekharl 
Il961h . thus possibly leading to longer cloud lifetimes. 
Two-di mensional numeric al models seem to support thi s 
notion (jKonz et all 120021 ; see also iPalotti et alJl2008f ). 
however, the dominant instabilities allow interchange 
modes in three dimensions, possibly with growthrates 
abov e the hydrodynamical ones (e.g. IStone fc Gardiner] 
120071 for the Ray leigh- Taylor instability). Currently 
available three-di mensional models o f HVCs including 
magnetic fields (Sant illan et alj I2004T ) are limited to 
clouds starting at z = 2.5 kpc above the plane and 
thus cannot address the question of cloud disruption over 
longer distances. Magnetic field estimates in HVCs are 
scarce - an estima te of 11.4 ± 2.4/iG in a s ingle HVC 
(iKazes et al.l fl991) could not be confirmed (jVerschuuri 
[HH: -0.1 ± 1.8/LtG for the same object). Indirect esti- 
mates based on cosmic ray confi nement quote a f ew /zG 
(|Valledll99d for the Galaxy, and lBeck et~aT1ll994l for th e 
starburst galaxy NGC 253; see also lGaensler et aLll2008f h 
With typical halo gas parameters, this would result in a 
plasma [3 of order unity, i.e. the fields might not be dom- 
inant. Clearly, the effects of magnetic fields need to be 
e xplored further. 

iVieser fc Henslerl (|2007bllah argue on the basis of high- 
resolution two-dimensional models of HVCs that heat 
conduction should play a substantial role in stabi- 
lizing the clouds against disruption by introducing a 
smooth transition layer between the hot and cold gas. 
Yet the effect seems to extend only the lifetime of a 
uniform-density cloud by a factor of two, while cen- 
trally condensed clouds (see Figs. 4a, 6a and 9a of 
IVieser fe Henslerl l2007bl ) are not strongly affected, pre- 
sumably because they already possess a smooth transi- 
tion layer, intrinsically weakening the dynamical insta- 
bilities. 

Viscosity could also play a role in stab ilizing the cloud 
again st dynamical instabilities (e.g. IChandrasekharl 
lT96l[ ) by introducing a finite shear layer and thus re- 
ducing the relative velocities essential for a rapid (and 
unconditional) growth of the Kelvin-Hclmholtz instabil- 
ity. Since viscosity is included in our numerical scheme 
at (resolution-limited) Reynolds numbers smaller than 
those expected for the Galactic halo, its stabilizing effect 
is actually over-estimated in our models. 

Given the distance constraints on the HVCs discussed 
here, our model HVCs are not massive enough to be 
self-gravitating, nor a re they expect ed to be dominated 
by dark matter (e.g. I Wakkeri [200l . Thus, we are ne- 
glect ing the otherwise stabilizing effect of self-gravity 
(e.g. iMurrav et aDll993T) and of external potentials (e.g. 
iMori fc Burkertll2000D . 

4.3. Consequences for Feeding the Galactic Disc 



Taking our estimated disruption timescales and dis- 
tances, we find that clouds at 10 kpc or less above the 
disk with HI masses < 10 4 5 M are highly unlikely 
to make it to the disk in the form of neutral hydro- 
gen clouds. Because of the limited mass range sampled 
(3 x 10 2 < M HI < 4.4 x 10 4 M Q ), we can only extrap- 
olate our model data to predict the evolution of larger 
HVCs such as Complex C (at « 10 6 M ; see Figure 6). 
These large complexes are often seen to consist of nu- 
mero us smaller clouds when observed at higher re solu- 
tion (|Schwarz fc Wakker|[200l iPutman et all 120021: Hsu 
et al. 2009), so much of their mass may also have diffi- 
culty making it to the disk in cold form. 

The lack of surviving clouds has implications for 
the accretion of gas from satellites. The Magellanic 
Stream is an example of a larger complex breaking 
up in to smaller clouds at the tip ( Stanimi rovic et al.l 
l2008f ) and via the numerous sm all clouds found along 
its length (jPutman et al.ll200 3b). Given the distance of 
the Magellanic Clouds (~ 60 kpc) the majority of the 
Stream is unlikely to mak e it to the disk in the form 
of cold HI clouds (see also iBland-Hawthorn et al.1 120071 
iBland- Hawthorn! l2008fl . Exclud ing the Magellanic Sys- 
tem. iGrcevich fc Putmanl (|2009f ) find that all satellites 
within 270 kpc of the Milky Way or Andromeda are 
devoid of HI gas. This gas was most likely primarily 
stripped at the satellite galaxy's perigalacticon, but this 
still means the gas was stripped from satellites at typ- 
ical distances of 20-120 kpc. Given these distances the 
majority of the gas stripped from satellites is unlikely to 
make it to the galaxy in cold form. If the fate of satel- 
lite galaxy gas is to ultimately provide the Galaxy with 
star formation fuel, it must first be integrated into the 
previously existing warm/hot halo medium. 

The destruction of cold clouds in the halo results in 
a complex multi-phase halo of warm gas intermingled 
with the hot and cold g ciS, clS also evident through 
the absorption line ob servations (|Sembach et alj 120031 : 
iBen Bekhti et al.l 120081) and pulsar dispersion measure- 
ments |Gaensler et al.1 [20081 although given the pulsar 
distances, most of the material they see should be part of 
a thick disk component. See also £|4.2.1[) . Even if a large 
fraction of the clouds do not reach the disk in the form of 
neutral hydrogen, they still are over-dense with respect 
to the backgr ound medium, and thus sink towards the 
plane (see also Bland-H awthorn et al.ll2007[) . One might 
speculate whether such remnants contribute to the pop- 
ulation of warm ionized medium (WIM) clouds in th e 
hot coronal gas (jRevnoldsl [l99l iGaensler et"aTll2008f) . 
Once they are within 1-2 kpc of the Galactic disk, the 
gas flows driven by feedback from m assive star form ing 
regions (Galactic fountain, see e.g. iBregmanl |2004| for 
a discussion) will most likely lead to compressions and 
recooling, eventually integrating the former HVC mate- 
rial into the disk as new star formation fuel. This could 
be easily tested by placing tracer particles in models of 
stratified Galactic disks including s upernova feedback, 
such as bv lJoung fc Mac Low! (|2006f ). 

The WIM phase is difficult to trace in our simulations, 
as it moves at slower velocities than the HI clouds due 
to being more diffuse, and thus leaves the simulation do- 
main rather quickly (within 39 Myr at 50 km s _1 ). Yet 
in models Hclbl3a and Hclbl5b (Fig. [4]) we were able 
to follow the cloud evolution long enough to observe its 
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nearly complete deceleration close to the plane, and the 
subsequent recooling of stripped material. These mod- 
els suggest that warm ionized gas stripped from a HVC 
can recool at low z (even without gas dynamics within 
the disk, as mentioned above), and thus re-form cold 
HI clouds. Because of their velocities, they would be 
classified as intermediate or low-velocity clouds (IVCs, 
LVCs) . These LVCs may be similar in nature to the dis- 
crete clouds recentl y detected in multiple directions by 
several HI surveys (|Lockmanl l2002t iStanimirovic et all 
l2006t iFord et al.ll2008[ ). These observed clouds are small 
and cold and are most likely embedded in the WIM ex- 
tending ~ 3 kpc above the plane. If they are a mixture 
of accreting gas from the halo and fountain material, the 
LVCs should have sub-solar metallicities and may them- 
selves show head-tail cloud structures as they sink into 
the Galactic disk. 

5. SUMMARY 

Mo tivated by the obse rved head-tail structure of HVCs 
(e.g. iBriins et alj|200d : Putman et al., in prep.; Fig. ^ 
and their potential to constrain the properties of the 
diffuse Galactic halo, we study numerically the disrup- 
tion of HVCs during their passage through the Galactic 
halo. Our experiments address the two possible extremes 
of HVC trajectories, namely (a) assuming a constant 
background halo density and pressure, and (b) assum- 
ing an isothermal, hydrostatic halo with the cloud falling 
straight towards the disk. These are our findings: 

Cloud Morphology and Evolution — In agreement with pre- 
vious numerical studies, the clouds develop a pronounced 
head-tail structure due to the interaction with the back- 
ground halo medium (Fig. [2]). Cooling does not neces- 
sarily stabilize the cloud, especially when the dynamical 
timescales are shorter than the cooling timescale (Fig. [3]). 
The velocity difference between head and tail is consis- 
tent with observed lags (Fig. 0|, and we suggest that ob- 
served radial cloud velocities should be correlated with 
the shape of the cloud and the head-tail lag. 

Cloud Disruption and Characteristic Distances — The lim- 
ited numerical resolution and our conservative parame- 
ter choice ( §4.2. 1|) renders the estimated cloud disrup- 
tion timescales and distances as upper limits. Therefore, 
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HVCs with HI masses < 10 4 5 M will lose their HI con- 
tent within 10 kpc or less (Fig.[6|). Extrapolating to more 
massive complexes (Fig. [7]), extreme ram pressures (e.g. 
7i/j > 10~ 2 cm -3 and v > 100 km s _1 ) or the fact that 
these complexes have abundant substructure (see §4.2. 1|) 
could lead to much of their HI mass being lost before 
reaching the disk. Material from the Magellanic Stream 
and other dwarf galaxies is not expected to survive the 
trip to the disk in the form of HI clouds. 

The Lower Halo and Upper Disk: Fueling Star Formation — 
Warm ionized material stripped from the HI clouds is still 
falling towards the disk (albeit at lower velocities). This 
material might co nstitute some of the extended layer o f 
warm ionized gas (|Revnoldslll993l ; iGaensler et alj [2008). 
Some of this material can actually recool and thus re- 
form HI cloudlets close to the disk. Because of their 
low velocities, such clouds would be identified as IVCs 
or LVCs. We argue that the HVC remnants still could 
merge with the Galactic disk as WIM clouds or via re- 
cooling close to the disk, and eventually fuel the star 
formation in our Galaxy via a warm cosmic rain. 
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